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The scattering amplitude for the recently discovered exactly solvable shape in- 
variant potential, which is isospectral to the generalized Poschl- Taylor potential, is 
calculated explicitly by considering the asymptotic behavior of the X\ Jacobi excep- 
tional polynomials associated with this system. 
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In recent years the ideas of Supersymmetric quantum mechanics (SQM) and shape 
invariant potentials (SIP) have greatly enriched our understanding of the exactly solvable 
potentials [1]. The search for the exactly solvable potentials has been boosted greatly due 
to the recent discovery of exceptional orthogonal polynomials (EOP) (also known as X n 
Laguerre and X n Jacobi polynomials) {2} y|. Unlike the usual orthogonal polynomials, 
these EOPs start with degree n > 1 and still form a complete orthonormal set with respect 
to a positive definite innerproduct defined over a compact interval. This remarkable work 
lead Quesne [4| to the discovery of two new SIP whose solution is in terms of Xi Laguerre 
and X\ Jacobi polynomials. Subsequently, a third SIP was discovered whose solution is 
also in terms of X\ Jacobi polynomials [5]. Subsequently, Odake and Sasaki constructed 
infinite sets of new SIP corresponding to all these three cases whose eigenf unctions are in 
terms of X n Laguerre and X n Jacobi polynomials jg]. 

It is worth pointing out that unlike the usual SIP, the newly discovered SIP are ex- 
plicitly h dependent. Further, all of them are isospectral to the well known SIP. Besides, 
while two out of the three newly discovered SIP have pure bound state spectrum, the 
third SIP which is isospectral to the generalized Poschl- Teller (GPT) potential, has both 
discrete and continuous spectrum. To the best of our knowledge, while the bound state 
energy eigenvalues and eigenfunctions have already been obtained in all the three cases, 
the scattering amplitude has still not been obtained in the case of SIP which are isospec- 
tral to GPT. The purpose of this note is to fill up this gap partially. In particular, in 
the present work we obtain the scattering amplitude for the newly discovered SIP whose 
solution is in terms of Xi Jacobi polynomial. 

The superpotential corresponding to GPT (which is on the half line < r < 00) is 
given by 

W GPT = A cothr - Scosechr B > A + 1 > 1 . (1) 

The potential VcpTij) = W^gpt^) - Wq PT {t) which follows from the above super potential 
is given by (h = 2m = 1) 

V GPT (r) = A 2 + [B 2 + A(A + l)]cosechr - B(2A + l)cosechr cothr . (2) 

whose bound state energy eigenvalues and eigenfunctions are well known. Remarkably, if 
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we consider 

2£?sinhr 2£?sinhr . . 

W = Wrpr-\ , (3) 

GPT 2B coshr -2A-1 2B coshr - 2A + 1 ' 1 ; 

then we find that even though the potential W 2 (r) — W'(r) is very different from Vqpt 

and given by 

x r / \ -r r 2(2A+1) 2[4B 2 - (2A + 1) 2 ] 

V(r) — VrPT-\ (4) 

yj GPT ^ 2B coshr - 2A - 1 (2Pcoshr - 2A - l) 2 ' K) 

the bound state spectrum is still the same and given by 

E v = A 2 - {A - vf , u = 0,l,..., 1 / max , (5) 

where A — 1 < v max < A. However, the eigenfunctions are now different and they are now 
given by 

AT (coshr - l)l^)(coshr + l)-*(*M> ( 
iMO = ^ 2gcoshr _ 2A _ 1 J (coshr) (6) 

where a = B - A - \ , /3 = -S - A—\. Here 

v _ 9 a +2 H / H(2A-2z,)(i? + A-z,+ i)r(i? + A-^-i) \ * 

V(s-i4 + v + i)r(s-i4 + i/-|)r(2>i-i/ + i); ' 1 j 

is the normalization constant, and P„+P is [y + 1) th-degree Xl- Jacobi Polynomial. 



The X\ Jacobi polynomial is related to the usual Jacobi polynomial as 



p(^) (r) = _± (a . _ 6)P^f ( r ) + u± r 1 V v " 2 I ; (8) 



where b = 

p— a 



Using OH]) the X\ Jacobi Polynomial -Pj+f "* (cosh r) can be written as 
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Pj+f (coshr) = — — {(b - coshr) (a + + 2v) + 2b] P^ (coshr) - 2Pj;lf (cosh 

2 (a + p + 2z/) L 

(9) 

Usual Jacobi polynomial p( a '^ (coshr) further can be written in terms of Hypergeometric 
function as : 

i„ R \, \ T(v + a + 1) „ 1 — coshr, 

P^ cosh r) = \ } F(u + a + P + l,-v,l + a; ) . (10) 

v\Y (1 + a) I 
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To get the scattering states for this system two modifications of the bound state wave- 
functions have to be made [3]: (i) The second solution of the Schrodinger equation must 
be retained - it has been discarded for bound state problems since it diverged asymptot- 
ically, (ii) Instead of the parameter v labeling the number of nodes, one must use the 
wavenumber k so that we get the asymptotic behavior in terms of e ±thr as r — > oo. After 
considering the second solution, Eq. (TlO|) becomes 



Pj^cosh; 



T(u + a + 1) 



+ C 2 



v\T(l + a) 
, 1 — cosh r 



dF{u + a + f3 + 1, -v, 1 + a 



1 — cosh r , 



F(y + (3 + 1, — v — a, 1 — a; 



1 — cosh r , 



2 v ' ' '2 

where Ci,C2 are arbitrary constants. Considering the boundary condition, i.e as r — > , 
( 1 ~ c ° !ihr ) _ y 0,^v(r) tending to finite, the allowed solution is 

Pi^cosh r) = T ±^±A Cx F{v + a + fi + 1, -„, 1 + a; 1^1) . (i 2 ) 
z/!l (1 + a) 2 

Now replace v by A + ik and use a + (3 = —2A — 1 We get 



P/?',%(coshr) = C 1 7 - r( ^,^ + 1 / 2) , , F(-A+ifc,-A-ifc,£-A+l/2; 



1 — cosh r , 



(A + ik)\V(B - A + l/2)' 



(13) 



Using Eq. PD in © we get (cosh r) = P ( ( ^ +1) (coshr). 

Thus the scattering state wavefunctions becomes 



Mr) = N k 



(coshr - l)K B ^)(coshr + 1)-*C*M) 

^A+ik+l l cosn r ) 



2B cosh T-2A-1 



(14) 



Now using the properties of hypergeometric function [8j i.e 

r(V)r(/3 - a) 1 



r(/3)r( 7 -a 
_ r( 7 )r(a-/3) 



P(/3,7-a,/3-a + 1: 



1 - z' 
1 



(15) 



and taking the limit r — >■ oo the fourth term of the hypergeometric equation vanishes. 
Finally we get the asymptotic form of (JHJ), which is given as 



lim if)k(r) — N} 



Ci2 



-2ik-3A 



[(2ik - l)aP + Qc] 



AB(2ik - 1) 



bP(l-2ik)2-* ik £ikr _ e _ ikr 



aP(2ik - 1) + Qc 



(16) 
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where 



T(B + ik + 1/2) T(B + ik- 1/2) 

— I A i ■i\lri/n A i 1 /o\ ' t> 



(A + ik)\T(B- A + 1/2)' ^ (A + life- l)!r(S-A + 1/2) : 



r(S - A + l/2)r(-2ifc) , T(B - A + l/2)T(2ik) 
a ■ 



T(-A - ik)T(B - ik + 1/2) ' + ik)T(B + ik + 1/2) 



r(B- A + l/2)r(-2zA; + 2) r(B - A + l/2)r(2iJfc - 2) 



r(-A - ik + 1)T(B - ik + 3/2) ' r(-A + ifc - l)r(S + ik - 1/2) : 



The asymptotic behavior for the radial wavefunction (for 1=0) is given by 
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\im Mr) - jj^ S i=^ ikr ~ e~' kr ] (17) 



From ( TT6l) and ( TT71) we get 

bP(l - 2ik)2~ 4ik 

Si=o 



aP(2ik -1)+Qc 

Using P, Q, a, b and c, we get after simple calculation (using F(n + 1) = nT(n)) 
_ qGPT [B 2 - (ik - 1/2) 2 ] 

^=0-^=0 [ jB 2_ ( , A;+1/2) 2 ] 

r(2ijb)r(-v4 - zA;)r(5 - zA; + l/2)2~ 4ifc [B 2 - [ik - 1/2) 2 ] 



(19) 



T(-A + ik)T(-2ik)T(B + ik + 1/2) [B 2 - (ik + 1/2) 2 ] 

Thus we notice that the scattering amplitudes for the two potentials (i.e. GPT and new 
Xi SIP) are different even though the bound state spectrum is identical for them. It will 
be interesting to see how the scattering amplitude change as we go from the potential 
here to the one whose eigenfunctions are in terms of X n Jacobi polynomials. 
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